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1.  INTRODUCTION 
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This  Final  Report  R930007-F  summarizes  activities  sponsored  by  the 
U.S.  A ray  Research  Office  under  Contract  No.  DAAG29-81-C-0033.  The  technical 
period  extended  from  15  August  1981  to  30  September  1982. 

Three  distinct  levels  of  study  were  involved  and  used  in  ARO-sponsored 
programs.  The  most  fundamental  level  addressed  quantum  transport  effects  in 
nanometer  scale  devices,  through  solutions  to  the  Wigner-Boltzmann  transport 
equation.  The  second  level  addressed  non-equilibrium  transport  on  a  submicron 
scale  through  solutions  to  the  Boltzmann-Block  transport  equation.  Here  such 
phenomena  as  velocity  overshoot  enter.  The  third  level  examined  transport  on 
a  near  micron  scale  through  solutions  to  the  semiconductor  drift  and  diffusion 
equations.  This  last  level  of  study  has  been  the  mainstay  of  semiconductor 
devices  for  the  last  thirty  to  forty  years. 

During  the  course  of  the  ARO  program,  several  papers  appeared  in  print  in 
which  quantum  transport  effects  were  studied.  The  motivation  for  the  study 
was  high  USER  (ultra  submicron  electronics  research)  interest  in  determining 
the  range  of  validity  of  the  Boltzmann  transport  equation  for  submicron  devices. 
Difficulties  are  suggested  by  the  fact  that  within  the  effective  mass  approxi¬ 
mation  the  thermal  deBroglle  wavelength  is  of  the  order  of  ultrasubmlcron 
dimensions  (e.g.,  250*A  for  GaAs).  Vith  this  as  motivation,  a  quantum  trans¬ 
port  equation  suitable  for  numerical  device  computation  was  developed.  Here 
the  equation  of  motion  for  the  Wigner  distribution  function  was  generalized 
to  include  scattering  contributions.  Balance  equations  for  density,  momentum 
and  energy  and  their  quantum  contributions  have  been  derived.  The  study  has 
demonstrated  that  quantum  corrections  are  non-negligible  when  the  carrier 
transit  lengths  are  of  the  order  of  the  thermal  deBroglle  wavelength.  Copies 
of  two  publications  on  this  subject  have  already  been  submitted  to  ARO. 

During  the  study  a  numerical  technique  was  developed  to  analyze  the 
small  signal  transient  characteristics  of  semiconductor  micron  length  and 
submicron  length  devices.  In  this  study  the  device  (either  two  or  three- 
terminal)  was  brought  to  steady  state  and  then  subjected  to  a  small  voltage 
pulse.  The  resulting  time-dependent  response  was  subsequently  determined 
self-consistantly,  and  then  analyzed  in  terms  of  frequency  dependent  Fourier 
components.  For  submicron  length  devices,  the  transient  analysis  was 
coupled  to  solutions  to  the  Boltzmann  transport  equation  and  the  small 
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signal  upper  frequency  limit  of  gallium  arsenide  two-terminal  devices,  operating 
under  transient  overshoot  conditions  was  determined.  The  results  were  incorporated 
into  a  major  review  article  and  were  published. 

The  small  signal  analysis  was  also  applied  to  the  semiconductor  drift  and 
diffusion  equation  to  obtain  the  small  signal  microwave  parameters  of  gallium 
arsenide  field  effect  transistors.  The  frequency  and  bias  dependence  of  these 
parameters  are  thought  to  directly  reflect  the  distribution  of  electrons  within 
the  device.  In  this  study  the  device  was  brought  to  a  steady  state  voltage 
level  and  then  subjected  to  a  small  voltage  pulse  on  the  drain  contact.  The 
resulting  time-dependent  device  response  was  obtained  through  two-dimensional 
solutions  of  the  relevant  semiconductor  equations.  The  time-dependent  response 
was  analyzed  into  frequency  dependent  (Fourier)  components  and  expressed  in 
suitable  microwave  ("Y")  parameter  form.  The  calculations  were  repeated  for  a 
pulse  on  the  gate  contact.  The  microwave  parameter  calculation  which  showed 
general  agreement  with  experiments  was  extremely  sensitive  to  the  electron 
distribution  within  the  device.  One  calculation  demonstrated  that  microwave 
parametric  measurements  could  reveal  the  presence  of  local  dipole  layers  of 
charge  within  the  device.  The  results  of  this  study  are  new  and  are  summarized 
in  this  report.  The  study  is  being  prepared  for  publication. 

Another  area  of  study  involving  the  drift  and  diffusion  equation  concerns 
calculating  the  noise  in  semiconductor  devices.  Particularly  in  the  case  of 
gallium  arsenide  devices,  there  is  a  large  body  of  experimental  evidence  showing 
an  increase  in  noise  prior  to  the  onset  of  Gunn  oscillations.  Indeed,  this 
result  led  some  early  investigators  to  speculate  that  high  field  cathode-to- 
anode  transit  time  oscillations  were  triggered  by  shot  noise.  In  contrast 
to  this,  theoretical  studies  show  a  decrease  in  FET  channel  noise  when  the  field 
is  in  the  region  of  negative  differential  mobility.  These  ostensibly  conflicting 
results  need  reconciliation.  One  phase  of  the  ARO  study  Involves  an  analytic 
study  of  the  effects  of  contact  originated  or  operational  originated  (e.g.,  an 
FET  at  high  voltage  levies)  nonuniform  field  profiles  on  the  noise  properties 
of  negative  differential  mobility  devices.  In  particular  with  the  "fixed 
cathode  boundary"  field  model  developed  to  explain  cathode  originated  Gunn 
instabilities,  the  mean  square  noise  voltage  per  unit  band  width  has  been 
obtained  analytically  using  the  "impedence  field  technique".  The  computed 
results  demonstrate  that  as  the  threshold  for  current  instabilities  is 
approached,  the  noise  component  increases.  This  is  the  first  theoretical 
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demonstration  of  enhanced  noise  as  the  instability  threshold  -  rather  than  the 
electron  transfer  threshold  is  approached.  The  results  were  summarized  in  a 
keynote  paper  presented  at  the  Surface  and  Interface  Conference  in  Trieste, 
August  1982,  and  are  discussed  in  this  report. 

One  of  the  significant  conclusions  of  the  tvo-terninal  noise  study  is 
that  in  an  active  aedia  noise  disturbances  propagate  and  siaultaneoualy  are 
amplified.  Generally,  amplification  within  the  active  region  of  an  FET  has  been 
accounted  for  by  introducing  an  effective  electron  teaperature.  To  aoae  extent 
this  approach  is  incomplete  insofar  as  it  Ignores  the  physics  of  growth  through 
transit  of  the  carriers.  To  overcome  this  difficulty,  we  have  generalized  the 
standard  noise  theory  of  FETs  within  the  framework  of  the  analytical  gradual 
channel  approximation.  For  analytical  purposes,  the  essential  features  of 
growth  are  retained  within  this  approximation.  In  addition  to  the  analytical 
noise  theory,  the  small  signal  study  demonstrates  that  the  dc  space  charge 
profiles  will  significantly  affect  the  noise  characteristics  of  the  device. 

We  have  qualitatively  coupled  the  two-dimensional  "Y"  parameter  calculation 
with  the  analytical  noise  study.  This  is  also  summarized  below.  It  too  is 
being  prepared  for  publication. 

Host  device  simulations  involving  the  compound  semiconductors  assume  a 
uniform  ambient  temperature.  Typically,  however,  the  actual  device  is  subject 
to  significant  temperature  variations  which  can  sometimes  permanently  alter 
device  properties  and  often  degrade  device  operation.  In  order  to  introduce 
design  changes  to  minimize  these  deleterious  temperature  effects,  the  detailed 
temperature  profiles  must  be  determined.  Under  the  present  ARO  contract  the 
semiconductor  device  numerical  code  has  been  generalized  to  Include  the  heat 
flow  in  gallium  arsenide  devices.  Additionally,  current  flow  now  includes 
transport  by  holes  as  well  as  electrons.  The  results  of  this  study  are 
summarized  below. 

A  new  study  was  initiated  to  determine  the  electrical  characteristics  of 
a  space  charge  limited  injection  field  effect  transistor.  The  motivation 
for  this  study  was  the  superior  performance  of  submicron  scale  injection 
FETs.  The  study  discussed  below  utilized  the  drift  and  diffusion  equations 
and  is  regarded  as  a  benchmark  calculation  to  which  submicron  calculations 
may  be  compared.  The  results  are  briefly  summarized  below. 
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The  technical  discussion  is  contained  in  the  next  four  sections.  Section  2 
consists  of  a  one-dimensional  analytical  study  of  the  effects  of  cathode 
originated  nonuniforanitles  on  the  noiae  properties  of  transferred  electron 
devices.  Section  3  contains  the  discussion  of  the  n  user  leal  "Y”  parameter 
determination  and  th«  analytical  formulation  of  the  FET  noise  problem.  Section  4 
includes  the  tempera  cure  dependent  study,  while  Section  5  contains  the  discussion 
of  the  injection  FET. 

A  list  of  recenc  ARO  sponsored  papers  is  also  included. 


2.  Noise  Properties  of  Transferred  Electron  Devices 


2.1  Introduction 

When  a  signal  is  applied  to  a  transferred  electron  semiconductor 
such  as  gallium  arsenide  or  indium  phosphide,  space  charge  waves  grow 
as  they  propagate  downstream  from  the  disturbance.  This  growth  mechanism 
is  indiscriminate,  and  both  desired  as  well  as  undesired  disturbances 
(noise)  are  amplified.  In  the  study  summarized  in  this  report,  two  types 
of  noise  calculations  were  considered.  The  first  type  was  one  dimensional, 
for  which  analytical  calculations  were  performed.  In  this  calculation 
the  goal  was  to  determine  whether  noise  calculations  could  provide  Insight 
into  the  nature  of  the  cathode  contract  conditions,  and  how  they  would  be 
modified  by  nonuniform  field  profiles.  The  calculations  represent  the 
first  detailed  analytical  study  of  noise  under  nonuniform  space  charge 
conditions  appropriate  to  the  onset  of  transferred  electron  oscillations. 

The  second  type  of  noise  calculation  is  associated  with  noise  in  FET's. 

This  calculation  is  a  combination  of  both  analytical  and  numerical  studies. 
Both  calculations  are  designed  to  demonstrate  the  effects  of  nonuniform 
profiles  on  the  noise  properties  of  devices. 

2.2  One-Dimensional  Calculations 

The  calculational  approach  to  one-dimensional  noise  problems  is 
discussed  with  reference  to  Figure  2-1.  Figure  2-1  displays  a  one-dimensional 
device  of  length  L  and  carrier  density  NQ.  A  small  imposed  signal  at  the 
source  of  the  device  manifests  itself  as  a  change  in  voltage  across  the 
device,  6^(t).  Similarly,  a  small  fluctuation  in  current,  of  magnitude 
A(X^)  at  the  point  X  ■  X^,  will  propagate  toward  the  anode  and  manifest 
Itself  as  a  change  in  voltage  ^^(t).  Shockley,  et  al  (1966)  considered 
two  such  disturbances:  one  where  the  fluctuation  resulted  in  an  Increase 
in  local  charge  density,  the  second  in  which  the  fluctuation  resulted  in 
an  equal  decrease  in  local  charge  density.  Each  disturbance  results  in  an 
impedance  change  Z(X^)  and  Z^).  Of  relevance  for  Xj  •  X^  +  AX,  is  the 
differential  impedance  and  "the  Impedance  field  defined  as 


vz 


Z(xp  -zcxp 

AX 
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in  which  terns  the  nean  square  noise  voltage  per  unit  hand  width  is  defined 
[Shockley,  et  al  (1966)3 
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In  the  above  N  is  the  carrier  density,  D  the  diffusivity,  and  e  the 
magnitude  of  the  electron  charge.  In  the  analytical  calculations  D  is 
taken  as  constant. 

Since  the  noise  calculation  describe  departures  from  steady  state, 
the  relevant  equations  are:  (1)  The  equation  for  total  current 
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e(NV  -  DVN) 
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and  (2)  Poisson's  equation, 

V2f  -  ~  (H  -  N  ). 

C  o 

The  electric  field  in  these  equations  is 

F  -  -Vf 

and  the  potential  drop  'fQ  is 

t„  ■/  ut 

O 

In  the  analytical  noise  calculations  the  underlined  terms  are  ignored. 

In  the  calculations  that  follow,  two  types  of  structures  are  considered: 
(1)  space  charge  limited  diodes,  and  (2)  diodes  with  thermally  generated 
carriers.  In  both  cases  the  effects  of  boundaries  on  the  noise  properties 
are  emphasized.  In  these  calculations  it  is  convenient  to  work  with  a  set 
of  dimensionless  variables  identified  in  Table  1. 


2.2- 4 

2.2- 5 
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In  dimensionless  terms,  the  relevant  equations  are  (in  one  dimension) 
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and 


2.2-9 


2.2-10 


In  equation  2.2-8,  the  second  differential  equation  is  relevant  to  space 
charge  limited  (SCL)  diodes.  (Note:  the  left  hand  part  of  equation  2.2-10 
is  in  ordinary  units) . 

In  examining  noise,  small  signal  perturbations  from  the  steady  state 
are  considered,  thus 

j*>30+«j  2.2-11 


f  =>  f  o  (x)  +  6f 
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where 


v  =>v  +  6v  *  v  +  p8f 
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dv  2.2-14 
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for  simplicity  v  is  taken  to  be  represented  by  three  linear  pieces,  as  in 
Figure  2.2. 

Ihe  relevant  small  signal  current  equation  is 
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for  either  doped  or  SCL  diodes.  In  the  above,  the  tine  dependence  is 
taken  as  e*w*.  The  solution  to  equation  2.2-15  for  both  uniform  and 
nonunifom  fields  is  part  of  an  extensive  literature  [see,  e.g.,  Shaw,  et  al 
1979)].  By  way  of  introduction  we  carry  this  solution  for  uniform  fields 
to  the  point  where  the  impedance  field  is  obtained  [see  also  Thim  (1971)]. 

It  is  then  generalized  for  nonuniform  fields. 

For  uniform  fields,  a  disturbance  originating  at  the  point  x^,  with 
amplitude  A(xj)  results  in  a  perturbed  field 
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where 
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The  potential  drop  arising  from  this  disturbance  is 
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And  the  impedance  field 
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which  is  equation  (6)  of  Thim  (1971). 

For  nonuniform  fields  the  solutions  are  more  complex.  For  both 
doped  and  SCL  diodes  the  solution  to  equation  2.2-15  is 

x  /  x 
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Differences  in  Che  results  of  SCL  and  doped  diodes  arise  after  integration. 
For  doped  diodes 


f 


2 

g(C)dC 


iutCx^.x^)  +  log 


j  /v(x_)-l) 


x2 

t^xVx2)  +0\0 

X1 


2.2-22 


2.2-23 


For  SCL  diodes 
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Integrating,  we  obtain,  for  doped  diodes 


x 


and  for  SCL  diodes 


x 
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to  O 


4kTR 
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For  fixed  cathode  fields,  as  in  the  case  of  gallium  arsenide  devices, 
a(u)  “  o.  For  fixed  cathode  conduction  currents  as  in  the  case  of  some 
indium  phosphide  devices  a(u)  -  6j(x  -  o)/iu>.  These  latter  results  will  be 
included  at  a  later  time. 

We  now  examine  the  field  response  when  a  disturbance  originates 
somewhere  within  the  device.  In  particular:  for  a  disturbance 
originating  at  a  point  x^,  and  within  either  the  ohmic,  NDM  or  saturated 
drift  velocity  region,  the  potential  drop  due  to  this  disturbance  is 
required.  For  the  doped  diode 

fiXx^A.u)  - 

*1 

with  6f(x  >  x^,x,«i>)  given  by  equation  2.2-28  with  a(w)**0  and  &j(u)  replaced 
by  A(xj).  Integrating  the  above,  and  including  a  disturbance  originating  at 
the  point  x 2  >  x^  leads  to  an  impedance  field 

3z(Xi)  ^-Xtfx^A)  2.2-31 

3x  “  X 

which  has  the  same  form  as  that  of  equation  (2.2-19).  For  SCL  diodes 


/* 


&{  (x>x^,xw)dx 


2.2-30 


3z(x!>  i.e~1>l,t(xlA)  2.2-32 

3x  "  im 

[see  also  Thornber  (1973)]. 

With  the  above  result,  we  are  in  a  position  to  evaluate  the  noise 
within  either  a  doped  or  space  charge  limited  diode.  We  first  illustrate 
the  result  for  uniform  fields  and  doped  diodes. 

For  uniform  fields,  with  D  m  1,  and  n  m  1 


&n/2* 


2 

-X(l-x1)/vQ 

1-e 
X 


2.2-33 


.  \ 

\/-vi/i0  \  -Wi0  * 

pie  cos  0-1 J-  «e  sin© 

'  '  (v2+u2)l 
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2.2-35 


(See  also  Thim  (1971).  In  the  above,  6  Is  the  transit  angle 


o 


We  note  that  for  ohmic  conduction  y  »  1,  and  for  t  sufficiently  large 


(w2)  4kTR 
x  or  m  o 

A0/2W  2 

1+  u 


For  conduction  within  the  NDM  region  y  <  0  and  the  mean  square  noise 
voltage  per  unit  bandwidth  increases  as  e+^*^o.  For  nonuniform  fields 
originating  within  and  remaining  within  the  NDM  region 


<«z> 

a572»a 


4kTR  ,  j 
_ o  .A.  o 

2.  2  i  A 

y_+u 

O 


t(OA) 


2y 
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which  bears  a  clear  similarity  to  equation  2.2-34.  There  are  important 
differences  between  the  two.  In  the  uniform  field  case  the  noise  depends 
critically  on  length,  in  which  the  longer  the  device  the  greater  the  noise. 
For  the  case  represented  by  equation  2.2-37,  the  noise  is  dependent  on  the 
relative  values  of  the  mean  carrier  velocity  j  and  the  velocity  of  the 
entering  carriers.  This  occurs  through  the  relation 


-yt(o,A)  ,  VV(A>  2.2-38 

*  jo-v(0> 

Thus,  as  J  -►  v(0)  the  noise  increases.  The  context  in  which  this  occurs 
is  contained  in  Figure  2-3  which  shows  a  nonuniform  field  profile  for  a 
ten  micron  long  gallium  arsenide  element.  Here,  as  the  bias  is  increased 
that  portion  of  the  element  within  the  NDM  region  increases,  and  this 
feature  leads  to  enhanced  noise.  Thus,  from  this  point  of  view,  the  content 
of  equation  2.2-37  is  similar  to  the  length  dependence  of  equation  2.2-34. 

For  the  diode  whose  structure  is  represented  by  figure  2.3,  which  is 
common  to  that  of  a  diode  just  prior  to  that  of  a  Gunn  type  oscillation 
the  mean  squared  noise  voltage  per  unit  frequency  range  is 


12 


2.2-39 


<(0  gtluh  .  <‘il  ,  /i-A  \ 

ZS717  1-Ha2  An/2*  1-Ha2  ' 

Rather  than  deal  with  this  equation,  we  follow  Thin  (1971)  and  examine 
the  noise  figure  F  of  a  high  gain  amplifier 

F  ~  1  +  M  2.2-40 


where 


M 


Afl/2*  4kTRQ | Re  z(»)| 


2.2-41 


and  Re  *(w)  is  the  real  part  of  the  device  impedance.  Z(w)  is  displayed 
in  figure  2-4  (see  Grubln  et  al  1975)  for  a  range  of  bias  values.  We 
note  that  z(ul  exhibits  zeroes,  and  that  the  separation  of  zeroes  is  an 
estimate  of  the  extent  to  which  the  nonlinear  element  sustains  electric 
field  values  within  the  NDM  region.  This  in  turn  introduces  structure 
into  F  [see  also  Sitch  et  al  (1976)J  which  is  shown  in  figure  2-5. 

The  noise  figure  thus  provides  two  pieces  of  information. 

First  the  noise  figure  Increases  with  increasing  bias,  reflecting  enhanced 
contributions  from  the  NDM  region.  This  increased  noise  prior  to  the  onset 
of  Gunn  oscillations  is  an  important  feature  of  Gunn  devices  (Gunn  1964). 
Secondly,  the  noise  figure  provides  a  means  of  tracing  the  extent  of  the 
NDM  region,  as  shown  in  figure  2-5.  As  the  bias  increases  the  edges  of 
maximum  F,  which  are  associated  with  ReZ  -*■  0,  move  closer,  reflecting  an 
increase  in  A.  This  narrowing  frequency  range  is  a  measure  of  transit 
across  the  NDM  region  where  gain  is  occurring.  The  significance  of  the  above 
unit  is  that  noise  measurements  may  be  utilized  to  provide  information  about 
the  structure  of  the  space  charge  layer  within  the  device.  Further  applica¬ 
tion  for  both  accumulation  layers  and  constant  cathode  conduction  current 
boundary' conditions  and  SCL  diodes  is  planned  for  the  future. 


Table  1 
Normalization 


Current  Density 

•>  -  «„•»,» 

Velocity 

V  -  V  v 
P 

Time 

T  -  T  t 
o 

Distance 

X  -  (V  T  )x 
P  o 

Diffusion 
D  -  (ykT/e) 


Carrier  Density 

N  -  N  n 
o 

Poetntial 

*  -  w* 

Field 

F  ”  V/pf 
P  ° 

Impedance 

Z  -  R(V  I  /L): 
o  p  o 

Angular  Frequency 
(2  -  T-1u 


Typical  Valves 

N  -  1015/cn3 
o 


Figure  2.1.  Schematic  of  small  signal  calculations.  A  small  current  disturbance 
at  source  results  in  potential  change  6i|i(T)  at  anode.  For  noise  calculations  a 
disturbance  at  x..  results  in  a  change  in  anode  potential  dijix. . 


ELECTRIC  FIELD  (f) 


Figure  2.2.  Three  linear-piece  representation  of  a  nonlinear  velocity-field 
curve.  In  normalized  variables  (see  Table  1),  the  velocity  as  well  as  the 
threshold  field  for  negative  differential  mobility  are  unity.  The  quantities  v> 
are  differential  mobilities,  expressed  in  normalized  units. 
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Figure  2.3.  Electric  field  versus  distance  profiles 
for  a  dimension  long  NDM  with  a  fixed  cathode  field 
Fc«  3.5Fj,  and  a  normalized  negative  differential 
mobility  of  -0.2.  Normalized  cathode  velocity  is  0.5. 
The  normalized  current  density  for  the  two  calculations 
are  0.485  and  0.488  respectively.  Note  that  the  NDM 
region,  represented  by  A  is  wider  at  the  higher  current 
density.  [From  Grubin  et  al  (1975)]. 
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Figure  2.4.  Small  signal  impedance  (in  multiples  of  Ro)for  a 
negative  differential  mobility  element  with  the  parameters  of 
figure  2.3.  The  abscissa  displays  both  angular  frequency  and 
transit  angle.  For  part  (a)  j*0.485,  for  (b)  jo*0.488. 

[From  Grubln  et  al  (1975)}. 
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NOISE  FIGURE  (dB) 


Figure  2.5.  Noise  figure  versus  transit  angle  keyed  to  normalized 
bias  current  j  -  J0/N0eV  ,  for  a  ten  micron  long  element  with  a 
depletion  layer  profile. 
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3.  Noise  in  Field  Effect  Transistors 


3.1  Introduction 

One  of  the  significant  conclusions  of  the  two- terminal  noise  study 
is  that  in  an  active  media  noise  disturbances  propagate  and  simultaneously 
are  amplified.  Generally,  amplification  within  the  active  region  has  been 
accounted  for  by  introducing  an  effective  electron  temperature  [Pucel,  et  al 
(1975),  Baechtold  (1972)].  To  some  extent  this  approach  is  incomplete  Insofar  as 
it  ignores  the  physics  of  growth  through  transit  of  the  carriers.  To  overcome 
this  difficulty,  we  have  generalized  the  standard  noise  theory  within  the 
framework  of  the  gradual  channel  approximation.  For  analytical  purposes, 
the  essential  features  of  growth  are  retained  within  this  approximation. 

This  is  discussed  below. 

In  addition  to  the  analytical  noise  theory,  it  is  clear  from  the 
previous  section  that  the  dc  space  charge  profiles  will  significantly 
affect  the  noise  characteristics  of  the  device.  Thus,  we  need  a  full 
set  of  bias  dependent  small  signal  parameters.  These  parameters  have 
been  calculated  numerically  and  are  discussed  in  Section  3.2.  Section  3.3 
formulates  the  noise  problem,  and  3.4  contains  the  analytical  contributions 
to  the  FET  noise  study. 

3.2  Calculation  of  the  Small  Signal  "Y”  Parameters 

This  calculation  is  a  numerical  one.  The  equations  solved  are  the 
two-dimensional  continuity  equation 

-e  +  div  •  J  -  0  3.2-1 

3T 

the  drift  and  diffusion  equation 

1  -  -e  J  N($,T)  $(?)  -  D(?)  VNO?»T)  )  3.2-2 

Poisson's  equation 

?2*  (S,T)  -  +  f 


^  (S,T)  -  No($)J 


3.2-3 
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and  the  relevant  circuit  equations.  In  the  above  $  denotes  the  two- 
dimensional  position  vector.  The  FET  configuration  is  shown  in  figure 
3-1. 

The  small  signal  calculations  are  obtained  as  small  time  dependent 
perturbations  of  time  Independent  steady  state  solutions.  Since  we  are 
treating  gallium  arsenide,  the  possibility  of  large  signal  domain  transit 
time  oscillations  must  be  considered.  Since  this  would  only  serve  to 
complicate  matters  and  render  the  concept  of  small  signal  parameters 
specious,  attention  was  concentrated  on  thin  "ND"  product  FETs  that  do 
not  sustain  instabilities  [Grubln,  et  al  (1980)]. 

The  steady  state,  or  dc  characteristics  of  the  device,  were  self- 
cons  is  tent  ly  computed  and  are  displayed  in  figure  3-2.  Typical  carrier 
density  profiles  are  represented  in  figures  3-3  and  3-4.  Figure  3-3 
displays  results  for  a  bias  sufficiently  high  to  generate  a  high  field 
domain  under  the  gate  contact.  Figure  3-4  shows  a  result  for  broad 
depletion  layer  under  the  gate  contact. 

Figures  3-3  and  3-4  and  their  respective  current  and  voltage  values 
represent  the  starting  point  of  the  calculation.  Subsequently,  a  small  signal 
square  wave  voltage  (see,  e.g.,  figure  3-5)  was  superimposed  on  the  drain  contact 
potential,  with  the  gate  contact  potential  fixed  at  its  steady  state  value. 

The  resulting  change  in  source,  gate  and  drain  currents  are  then  computed 
subject  to  the  constraint 

41  (T)  -  41,(1)  +  41  (T)  3-2-4 

sag 

The  gate  and  drain  current,  and  the  drain  potential  are  then  Fourier 
analyzed.  The  ratios  for  a  fixed  gate  potential  are  identified  as  the 
admittance  parameters: 


tI2(0)  -  4Ic(n)/4*D(Q) 

3.2-5 

Y22(Q)  -  4ID(n)/4*D(fl) 

3.2-6 

A  similar  exercise  is  performed  for  a  perturbation  on  the  gate  contact,  with 
the  identification  of  two  additional  admittance  parameters 
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Yn(fl)  -  6Ig(Q)/6*g(Q)  3.2-7 

T21(Q)  "  «ID<0>/«*C(0)  3'2"8 

The  "Y"  or  admittance  parameters  are  dependent  on  the  space  charge 
and  potential  profiles  and  are  thus  bias  dependent.  In  addition  it  is 
generally  assumed  that,  about  a  given  bias  point,  the  small  signal 
currents  add  in  a  linear  way;  i.e., 

6ic(n)  -  Yn(Q)  6*g(Q)  +  Y12(n)  «*D(a) 

6ID(Q)  -  Y21<n)  6*G(n)  +  Y22(n)  6*D(n) 

Calculations  of  the  "Y"  parameters  represent  the  core  of  the  small 
signal  calculation.  The  calculation  takes  on  added  significance  when 
placed  in  the  context  of  an  equivalent  circuit  model.  A  purely  formal 
manipulation  of  equations  3.2-9  and  3.2-10  provides  the  necessary  direction 

4IC  -  «11  +  *12>  ‘*c  -  Y12  <**C  -  5V  3-2-n 

8I0  ’  (Y21  *  Y12>  +  l12  <5*0  -  «V  +  (,22  +  Y12>  S*D  3’2’12 

with  the  equivalent  circuit  as  shown  in  figure  3-6. 

The  terms  Yll  +  Y12  are  generally  identified  as  source-gate  admit¬ 
tance  parameters, while  ~Y^2  is  referred  to  as  the  gate-drain  admittance. 

In  the  discussion  that  follows  we  have  attempted  to  construct  the 
simplest  type  of  equivalent  circuit  consistent  with  the  numerical  calculations. 
The  essential  features  of  the  calculations  are  represented  by  the  equivalent 
circuit  of  figure  3-7  where  represents  channel  resistance  in  the 
drain  loop  away  from  the  gate  contact.  The  equivalent  circuit  of  figure  3-7 
will  require  modifications  which  will  be  discussed  in  the  conclusions,  but 
it  represents  the  essential  features  of  the  calculation.  Identifying 
zeroth  order  "Y"  parameters  by  the  subscript  "o",  the  "Y"  parameters  are 
obtained  from  the  equations 


3.2- 9 

3.2- 10 
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3.2-13 


“*  '  I110  «*G  -  *8  4Ig>  +  ’no  <%  -  Sd  4I„> 

4I„  -  ’’210  <4*C  -  “g  4V  +  Y220  <4*d  -  *,  4Id> 


Defining 


and 


det  Y  2  Y110  Y220  -  Y120  Y21q 


D  "  1  +  Rg  Y110  +  Rd  Y220  +  Rg  Rd  det  Y 


the  "Y"  parameters  are 


’ll 


.  Y110  +  Rd  det  Y 


Y12  “-IP 


Y21  '  210 
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220  +  it 


R  det  Y 


For  our  immediate  discussion  it  is  sufficient  to  consider  the  case 
Rg  ■  0;  thus 

Y..  -  Yno  -  Y120  Y210  Rd 
1  +  Y220  Rd 
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120 


1  +  Y220  Rd 


‘21 


‘210 


1  +  Y220  Rd 
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‘220 


1  +  Y220  Rd 


3.2- 14 

3.2- 15 

3.3- 16 

3.2- 17 

3.2- 18 

3.2- 19 

3.2- 20 

when 

3.2- 21 

3.2- 22 

3.2- 23 

3.2- 24 
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With  reference  to  figure  3-7,  the  zeroth  order  "Y”  parameters  are 


Y120  +  Y220  “  GDS 

where  G^g  is  the  drain  conductance; 

Y210  “  Y120  “  Gm 


3.2-25 


3.2-26 


where  G  is  the  transconductance,  and 

YU0  +  Y120  ■  +  J“CSS  V  3-2-27 

Thus,  the  central  circuit  parameter  to  identify  is  the  gate-drain 
admittance  Y^q*  This  parameter  contains  the  essential  physics  of 
domain  formation.  At  low  drain  bias  levels,  or  rather  a  bias  level 
in  which  high  field  domains  are  not  present,  Y^q  is  adequately  repre¬ 
sented  by  the  capacitance  as  shown  in  figure  3-7.  At  high  bias 
levels  is  replaced  by  the  element  shown  in  figure  3-8,  where  the 

conductance  at  high  bias  levels  will  exhibit  a  small  region  of  negative 
differential  conductivity.  In  the  qualitative  discussion  that  follows, 
we  will  ignore  the  real  part  of  Y^q  insofar  as  it  modifies  the  equivalent 
circuit  representation,  although  its  contribution  is  incorporated  i.n  the 
numerical  contributions.  The  bias  dependence  of  Cgd  will,  howeVet f  enter 
prominently  in  the  following  discussion.  In  terms  of  figure  3-7 


Y110  "  -JnCgs  "  ^Cgd 
Y120  *  jnCgd 
Y210  "  Gm  +  J°Cgd 
Y220  "  Gds  -  ^Cgd 


3.2-28 


3.2-29 


3.2-30 


3.2-31 
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Thus 


Tu  -  Ad  (1  *  Rd  Gl> 


-  JO  K.  +  C»d  (1  *  Bd  Gl)  U  +  »j  gj,)|  3.2-32 


where 


G 


1 


ds 


+  G 


3.2-33 


and 


» ■ (i  +  «4.  Rd>2  +  °2  *5  cg2> 


3.2-34 


Also, 


T„  -  Art  Rd  3-  Cgd  (1  *  Rd  Cd.> 

*■2  n  TV 


3.2-35 


T  -  Gn  (1  +  Rd  W  "  °2  CRd  Rd  +  j0Cgd  (1  +_Rd_V  3.2-36 
21  -  D 


and 


t22  -  Cd.  <  1  +  *d  Cd»>  *  Ad  Rd  -  J°c,d 


3.2-37 


Selective  "Y"  parameter  versus  frequency  are  displayed  in  figures  3-9 
through  3-12  for  a  range  of  drain  bias  values  and  a  moderate  value  of  gate 
bias.  Figures  3-13  and  3-14  show  the  effects  of  increasing  the  gate  bias. 

In  these  figures,  the  admittance  parameters  are  expressed  units  of  Gq,  where 


N  epo  HW 
o 


3.2-38 


the  frequency  is  in  multiples  of  f 


0.542  GHZ 


3.2-39 


and  Che  bias  levels  are  in  multiples  of  F^L,  where  is  the  threshold 
field  for  negative  differential  mobility  in  gallium  arsenide.  Some 
broad  general  features  are  clear  from  equations  3.2-32  through  3.2-37. 

Consider  first.  Re  Y^.  It  is  approximately  zero  at  low  frequency 
values  and  shows  an  increase  in  value  with  increasing  frequency.  The 
imaginary  part  is  always  negative.  Y11  displays  Interesting  structure 
as  domains  form.  (We  note  that  the  space  charge  distribution  accompanies 
each  figure.)  Next  consider  Y^2.  1°  Che  absence  of  domain  formation  Re  Y^2 

is  small  under  dc  conditions  and  then  everywhere  negative.  Under  domain 
formation  there  is  a  significant  change  in  Re  Y^>  requiring  the  addition  of  a 
small  negative  conductance,  as  in  figure  3-8.  For  this  case 
Re  Y12  reads 


8 »  <i  +  Rd  Gds)  -  £12  C  2  Rd 


3.2-40 


where  gm  is  the  parallel  conductance.  Thus,  under  weak  domain  formation 
Re  Yj^  <  0.  Of  course,  there  is  a  frequency  dependance  to  gm  and  if  gm 
becomes  positive  at  high  frequencies.  Re  Y^  may  become  positive.  The 
possibility  also  exists  that  Gdg  will  become  negative  when  domains  are 
present  and  again  change  the  sign  of  Re  Y12.  The  quantities  Y21  and  Y22 
shows  the  general  trends  of  equations  3.2-36  and  3.2-34,  but  also  appear 
to  show  some  unusual  structure  when  high  field  domains  are  present. 

The  details  of  this  structure  are  discussed  below.  As  the  discussion 
progresses,  necessary  changes  to  the  equivalent  circuit  will  become 
apparent . 

We  begin  with  the  detailed  discussion  of  the  "Y"  parameter  with  an 
examination  of  the  drain  components  Y^2  and  Y22.  An  examination  of 
these  Y-parameters  for  a  fixed  gate  bias  of  -0.1  normalized  units 
(figures  3-9  through  3-13)  indicates  that  at  low  values  of  drain  bias, 

|lm  Y22|  <  | 1m  Y12| .  This  situation  remains  until  a  high  field  domain 
forms  within  the  channel  and  |lm  Y22|  >  |  Im  Y^2|.  The  situation  in  which 
the  drain  bias  is  fixed  and  the  gate  bias  is  increased,  result  in  movement 
of  the  depletion  layer  toward  the  bottom  of  the  channel.  It  is  seen  that 
|lm  Y22|  >  |lm  Y12| .  Thus,  the  drain  element  YQ  appears  to  require  the 


25 


inclusion  of  a  capacitor.  In  any  case,  at  high  bias  levels, 

|lm  ^22 !  **  I ^in  Yljl  88  suggested  by  equations  3.2-35  and  3.2-36.  Further, 
the  frequency  dependence  and  bias  dependent  nature  of  both  Y^  and  Y^ 
suggest  a  rich  space  charge  dependency.  Consider  the  broad  frequency  depen¬ 
dence  of  Y^*  The  first  point  we  note  is  that  Re  at  low  bias  levels  is 
negative.  This  is  accounted  for  by  the  presence  of  the  parasitic  resistance 
associated  with  those  portions  of  the  semiconductor  element  not  directly 
under  the  gate. 

Apart  from  the  frequency  dependence  of  the  real  part  of  Y^*  t*ie 
most  dramatic  changes  are  those  associated  with  the  imaginary  path. 

As  seen  in  figure  3-15,  the  capacitance  of  the  system  undergoes  a 
percipitous  drop,  a  feature  also  seen  experimentally  [Englemann,  et  al  (1977)] 
figure  3-16.  The  broad  features  of  this  result  were  discussed  by 
Englemann,  et  al  (1977),  and  is  represented  in  figure  3-8.  Essentially, 
the  presence  of  the  high  field  domain  is  cutting  off  the  gate-drain 
coupling. 

He  recall  that  represents  the  gate-to-channel  capacitance  on 
the  drain  side  of  the  channel  as  indicated  in  figure  3-17,  and  is  a 
measure  of  the  change  in  channel  depletion  charge  as  a  result  of  changes 
in  drain  bias.  Now,  while  "*  capacitance  change  is  expected  as  the  edge 
of  the  depletion  layer  moves  toward  the  bottom  of  the  channel,  the 
enhanced  drop  in  capacitance  suggests  that  most  of  the  modulation 
is  across  the  dipole  layer.  This  is  illustrated  in  figure  3-15. 

There  are  basically  two  sets  of  data  in  figure  3-15.  The  bold 
line  is  primarily  for  a  moderate  value  of  gate  bias.  The  most  signifi¬ 
cant  drop  in  capacitance  occurs  when  a  domain  forms.  The  situation  when 
the  Increase  in  gate  to  drain  voltage  is  a  result  of  an  increased  gate 
bias,  does  not  result  in  the  dramatic  decrease  in  capacitance.  For  the 
latter  case,  dipole  domain  do  not  form.  He  note  that  the  qualitative 
features  of  this  calculation  are  in  agreement  with  the  experimental  results 
of  Englemann,  et  al  (1977) . 

For  the  parameter  Y221  we  simply  note  that  at  low  bias  levels 
Gds  >  0,  whereas  at  high  values  of  bias  the  results  strongly  suggest  a 
frequency  dependent  region  of  small  signal  negative  conductance. 


We  next  consider  and  note  again  that  to  account  for  Che  high 
frequency  high  bias  behavior  of  It  is  necessary  that  G^s  exhibit  a 
frequency  dependence  which  allows  for  a  range  of  snail  signal  negative 
conductance.  With  regard  to  In  Y,  ,f ,  which  we  write  as 

I®  YH  "  cu  <G)  3.2-41 


we  note  from  figures  3-18  and  3-15  that  at  low  bias  levels  is 
somewhat  greater  than  twice  C^.  At  high  bias  levels  when  domains  form 
Cn(&)  is  at  least  an  order  of  magnitude  greater  than  C^.  Under  low 
or  moderate  gate  bias  levels  C^(Q)  does  not  exhibit  a  precipitous  drop 
in  value.  Rather,  at  first  the  capacitance  decreases  corresponding  to  a 
movement  of  the  depletion  layer  toward  the  bottom  of  the  channel.  This 
was  also  the  initial  behavior  of  the  gate  to  drain  capacitance.  Further 
increases  in  drain  bias  result  in  domain  formation  and  space  charge 
injection  into  the  depleted  zone.  The  effective  capacitance  shows  a  cor¬ 
responding  increase.  This  is  displayed  in  figure  3-18.  For  the  situation 
where  the  net  voltage  increase  is  due  to  an  increase  in  gate  bias,  where 
no  domains  form  there  is  the  expected  drop  in  capacitance,  as  reflected 
in  the  gate  to  drain  capacitance  without  domains.  This  is  also  shown 
in  figure  3-18.  The  experimental  situation  shows  broad  agreement  with 
the  numerical  results  and  is  displayed  in  figure  3-19. 

The  two  remaining  items  of  interest  here  are  theoretical/experimental 
comparison  of  the  transconductance,  and  current-gain  cutoff  frequency. 
Figure  3-20  displays  the  transconductance.  Re  Y21,  versus  VDS  -  Vgg,  for 
two  values  of  gate  bias.  We  see  the  presence  of  near  saturation  in  both 
sets  of  data.  This  is  also  seen  experimentally.  Our  data  does  not  extend 
to  high  enough  drain  bias  levels  to  determine  whether  a  corresponding 
decrease  in  transconductance  occurs.  The  experimental  results  are  shown 
in  figure  3-21. 

The  current-gain  cutoff  frequency  is  obtained  from  the  expression 
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and  Is  shown  in  figure  3-22.  It  is  seen  that  saturation  in  f^,  occurs 
under  the  presence  of  domain  formation.  The  decrease  in  f^  at  the 
higher  drain  bias  levels  appears  to  be  associated  with  an  increase  in 
the  source-gate  capacitance.  Experimental  observations  are  displayed 
in  figure  3-23,  and  again  show  qualitative  agreement  with  theory. 

The  results  of  the  above  study  show  a  clear  indication  of  the  role 
of  the  high  field  domain  on  the  small  signal  characteristics  of  the  FET. 
These  will  be  used  in  the  qualitative  discussion  of  noise  contributions 
to  the  FET  in  the  next  section.  One  of  the  more  interesting  results  is 
the  presense  of  negative  values  of  Re  Y^.  This  requires  further 
examination. 


3.3  Circuit  Representation  of  Kolse  Sources 

As  discussed  in  section  3.2,  the  small  signal  properties  of  the 
FET  can  be  expressed  in  terms  of  a  set  of  "Y"  parameters,  where  when 
noise  is  included  the  circuit  is  as  shown  in  figure  3-24.  For  a 
network  of  this  type  it  has  been  established  that  the  noise  character¬ 
istics  of  the  network  can  be  completely  specified  in  terms  of  two  noise 
generators  [see  e.g.,  Talpey  (1959)]  as  in  figure  3-25  or  equivalently 
figure  3-26.  In  figure  3-26 


6!  -  6It  -  Yu  «ld/Y21 


3.3-1 


« *  -  «y*21 


3.3-2 


As  discussed  earlier,  a  quantitative  measure  by  which  the  FET  noise 
characteristics  are  generally  judged  is  the  noise  figure.  The  noise  figure 
is  defined  as  the  ratio  of  the  output  noise  power  of  the  actual  device  to 
the  output  noise  power  that  would  be  obtained  from  a  hypothetical  device, 
identical  in  all  respects  except  that  it  contains  no  noise  sources.  The 
output  noise  in  the  ideal  case  is  then  amplified  thermal  noise  from  the 
signal  source  connected  to  the  input  terminals.  Thus 
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3.4  Calculation  of  Noise  Sources 


"  3.4.1  Introduction 

The  quantities  of  interest  in  evaluating  the  noise  figure  are 
the  noise  current  sources  and  the  circuit  "Y"  parameters.  The  "Y" 
parameters  are  computed  self-consistently.  The  noise  parameters  are 
not.  Instead,  we  have  sought  to  develop  conceptual  ideas  associated 
with  the  propagation  of  noise  sources,  which  we  have  seen  tend  to 
dominate  two-terminal  device  behavior.  Thus,  while  the  frequency  and 
bias  dependence  of  the  noise  current  generators  are  dependent  on  the 
circuit  elements  and  should  be  obtained  self-consistently,  we  will  rely 
heavily  on  a  small  signal  formulation  that  is  consistent  with  the  gradual 
channel  approximation.  It  may  be  argued  that  the  gradual  channel 
approximation  will  not  properly  describe  the  electrical  behavior  of  gal¬ 
lium  arsenide  field  effect  transistors,  and  this  author  would  agree. 
Nevertheless,  from  the  view  point  of  highlighting  important  physical 
concepts  it  remains  extremely  valuable.  It  is  from  this  latter  point 
of  view  that  it  is  used  below. 

3.4.2  Critique  of  Noise  Calculations,  Nonlocal  Equations 

With  reference  to  figure  3-28,  the  total  channel  current  between 
source  and  drain,  ignoring  diffusion  contribution  is 

Id  -  (eNV  +  e  ||)  W(H-h(x))  3.4-1 

In  the  absence  of  any  gate  current  1^  is  independent  of  x,  although  in 
general  1^  is  position  dependent. 

The  approach  of  this  author  to  performing  noise  calculations 
was  expressed  in  the  one-dimensional  section.  It  was  assumed  that  a 
local  fluctuation  somewhere  within  the  device  caused  a  change  in  current, 
and  that  the  resulting  time-dependent  behavior  was  described  by  a  small 
signal  perturbation  equation.  For  the  FET  study  we  assume,  therefore, 
that  prior  to  the  fluctuation  there  is  no  gate  current,  while  after  the 
fluctuation,  gate  current  can  flow.  Within  the  spirit  of  the  small 
signal  analysis  we  take 
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h<x,T)  ->  h  (x)  +  6b 
o 

H(x,T)  ->  H  (x)  +  6N 
o 

V  ->  V  +  u6F 

O 

F(x,T)  ->  F  (x)  +  6F 
o 
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where 


o 


Incorporating  Poisson's  equation  and  rearranging  equation  3.4-1  we  obtain 
the  small  signal  equation: 


do 


+  61 , 


-Id 


W(H-ho(x)-6h(x)) 


W(H-ho(x) 


i-  ♦ 

\3T  o  JX  +  t  ] 


6F  3.4-4 


where 


do 


W(H-ho(x)) 


N  eV  (x) 
o  o  7 


3.4-5 


To  place  equation  3.4-4  in  a  form  useful  for  analytical  work,  we 
introduce  the  depletion  layer  approximation  for  h(x).  If  represents 
the  potential  along  the  channel,  and  6^  is  the  potential  on  the 
gate  contact 
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1^(1)  is  directly  lntegrable: 

3/2  i3/2 
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1^(2)  requires  numerical  integration.  The  significance  of  this  separa¬ 
tion  is  that  1^(2)  is  zero  when  NQ  is  spatially  Independent;  an  assumption 
made  in  virtually  all  analytical  FET  calculations  [see  e.g.  Pucel,  et  al 
(1975)].  Finally,  6h(x)  is  given  by 


6h(x) 
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Writing  the  small  signal  equations  in  terms  of  the  reduced  potential  $, 
we  obtain 
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In  the  above 
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Simllarlly  for  Sid.  We  have  also  used  the  Table  1  normalization  for  t 
and  n(x). 

Equation  3.4-17  is  the  nonlocal  equation  governing  small  signal 
transport  within  an  FET.  It  is  more  general  than  hitherto  proposed, 

[see  e.q.  Pucel,  et  al  (1975),  Van  der  Zeil  (1963)  1  insofar  as  it  contains 
the  seeds  of  a  propagating  disturbance.  For  example,  if  the  underlined 
term  on  the  left  hand  side  of  equation  3.4-17  is  ignored,  the  equation  is 
essentially  that  of  Section  2  in  which  the  one-dimensional  transport 
equation  was  studied.  Here  it  is  seen  that  a  disturbance  within  the  FET 
will  grow  if  the  field  were  within  the  NDM  region,  while  decay  would  occur 
for  transport  within  the  ohmic  or  lossy  region.  The  standard  calculation. 
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however,  ignores  the  underlined  terms  on  the  right  hand  side,  assumes 
transport  is  local  within  the  ohmic  region  and  that  n(x)  -  1.  Under 
these  more  common  assumptions  the  governing  equation,  with 

‘do  ■  -L  (»  &  • 

is 

61d  *  - L  h  j6**1  -V*  -  *g)| 

See  e.q.  Van  der  Ziel  and  Ero  (1964),  equation  14. 

In  the  absence  of  numerical  calculations,  the  question  that  must 
be  asked  in  assessing  the  standard  approximations  to  equation  3.4-17  is 
how  significant  are  the  propagating  solutions  and  how  well  does  equation 

3.4- 20  represent  the  essential  physics.  These  questions  are  Independent 

of  whether  the  dc  gradual  channel  approximation  is  extended  into  the  region 
of  negative  differential  mobility. 

In  the  one-dimensional  noise  study  it  was  seen  that  disturbances 
within  the  ohmic  region  decay  very  rapidly  away  from  their  origin.  For 
this  case,  equation  3.4-20  is  likely  to  be  an  excellent  approximation 
for  a  zeroth  order  charge  density  that  is  approximately  uniform.  For 
disturbances  within  the  NDM  regime  the  full  perturbation,  equation  3.4-17, 
is  needed. 

3.4.3  Drain  Noise  Current;  Qualitative  Effects  of  Nonlocality 

We  now  consider  an  application  of  these  equations  to  the  problem  of 
determining  the  equivalent  noise  source  in  the  FET.  Consider  equation 

3.4- 20,  within  the  framework  of  the  impedance  field.  With  reference  to 

figure  2-1,  we  assume  a  small  fluctuation  in  current  is  inserted  at  x^ 
and  removed  at  Xj.  The  effect  of  the  current  pulse  at  x^  is  to  cause  a 
voltage  change  at  the  drain  of  magnitude  such  that 

"di(l  -  *i> '  5*dia  -  V*d  -  V  3-4-21 
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To  integrate  equation  3.4-29,  note  is  taken  of  the  fact  that 
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where 
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Here  g^(x,y)  and  f^(x,y)  represent  the  numerator  and  denominator, 
respectively  of  equation  3.4-32.  [See  e.g.  Van  der  Ziel  (1963),  equation  17)] 
where  P  is  of  the  order  of  magnitude  of  unity.  The  above  results  are  the  low 
frequency  results,  and  ignore  nonlocality.  The  interest  here  is,  however, 
in  the  propagating  terms.  The  discussion  below  is  concerned  with  this 
point  and  the  approach  is  qualitative.  We  begin  in  a  manner  that  is 
analogous  to  the  one-dimensional  noise  problem. 

For  the  one-dimensional  noise  problem,  and  with  uniform  fields, 
the  impedance  field,  in  dimensioned  units  is 


where 


R 

o 

1  - 

exp  -  [iQ 

vhQ 

L  -  X 

3Z 

T 

V 

ax 

Lt 

u/w. 

- 

m  + 

0 

T 

T 

E 

N  ey 

0  o 

R  C 
o  o 

3.4-33 


3.4-34 


For  low,  or  zero  frequency 
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For  disturbances  within  the  ohmic  region,  the  exponential  contribution  is 
negligible  and  with  y  ■  yQ 
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With  the  exception  of  a  variable  cross-sectional  area  this,  of  course, 
is  similar  to  the  results  of  equation  3.4-25.  The  situation  with  the 
drain  current  noise,  as  dictated  by  equation  3.4-17  is  that  the  impedance 
field  should  have  a  form  qualitatively  similar  to  that  of  equation  3.4-33 
and  that  for  the  FET,  equation  3.4-25  should  be  generalized  to  read 
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where  x  ^  is  a  generalized  "RC"  time  constant  to  account  for  the  fact 
that  the  cross-sectional  area  is  variable 

If  equation  3.4-37  is  intuitively  accepted  as  a  reasonable  repre¬ 
sentation  of  the  nonlocal  contributions  then  modifications  to  equation 
3.4-31  are  in  order.  To  introduce  these  modifications,  we  turn  for 
guidance  again  to  the  one-dimensional  case  concentrating  on  equation  2.2-33. 
This  equation,  under  dc  conditions,  yields  for  the  mean  square  current  per 
unit  bandwidth. 
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The  above  expression  is  for  a  disturbance  within  the  negative  differential 
mobility  region.  The  above  equation  indicates  that  a  fluctuation  within 
the  NDM  region  will  be  amplified,  as  expected.  Now  .within  an  FET  only  a 
fraction  of  the  device  may  be  expected  to  be  within  the  NDM  region.  Thus, 
qualitatively,  if  <  A  >  represent  an  average  length  of  the  NDM  region, 
we  expect  that  equation  3.4-31  will  be  modified  to  read 
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In  the  above,  the  length  dependance  of  A  is  determined  by  replacing  L 
in  equation  3.4-40  by  <  A  >.  The  meaning  of  equation  3.4-41  is  that  in 
the  absence  of  amplification  within  the  NDM  region  noise  is  given  by  the 
well  established  classical  expression.  In  the  presence  of  NDM  gain  over 
a  distance  <  A  >,  the  excess  noise  over  and  above  that  associated  with 

A 

classical  noise  over  a  distance  <  A  >,  P<  ^  >t  is  given  by 
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In  an  FET  and  within  a  region  of  negative  differential  mobility  the 
electron  velocity  can  change  by,  in  many  cases  a  factor  of  two  in  going 
from  one  region  to  another.  Thus,  it  is  not  clear  from  this  expression 
how  a  current  increase  will  enhance  the  noise,  as  seen  experimentally, 
if  all  interpretation  is  through  velocity  increases.  On  the  other  hand, 
numerical  calculations  reveal  that  the  length  <  A  >  of  the  NDM  region 
Increases  with  Increasing  drain  bias.  This  will  increase  A  and  the 
noise. 


3.4.4  Current  Noise;  Quantitative  Formulation  Including  Nonlocality 

In  many  ways  the  formulation  of  3.4.3  is  more  general  than 
the  approach  which  assumes  a  variable  electron  temperature  within  the 
channel  [Pucel,  et  al  (1975)].  Certainly  carriers  in  the  central  valley 
of  gallium  arsenide  experience  an  increase  in  electron  temperature  prior 
to  and  during  transfer.  But  once  in  the  satellite  valley,  the  mean 
carrier  temperature  decreases.  This  feature  is  reflected  in  the  dif¬ 
fusion  coefficient,  whose  field  dependence  is  representative  of  only 
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moderate  mean  temperature  variations.  The  crucial  component  of  the 
electron  transfer  mechanism  is  gain  through  negative  differential  mobility. 
This  result  is  expressed  qualitatively  in  equation  3.4-41.  It  is  worth¬ 
while  noting  that  if  the  mechanism  of  negative  differential  mobility  is 
ignored  and  all  enhanced  noise  is  due  to  variations  in  electron  temperature, 
then,  fundamentally,  silicon  should  be  qualitatively  similar  to  gallium 
arsenide.  In  silicon,  high  fields  also  result  in  increased  values  of 
electron  temperature,  but  there  is  no  gain  due  to  negative  differential 
mobility. 

Equation  3.4-37  indicates  the  presence  of  transit  time  effects 
associated  with  the  region  of  negative  differential  mobility.  In  view 
of  the  fact  that  this  region  is  often  shorter  than  such  critical  dimensions 
as  the  gate  to  drain  spacing  [Grubin,  et  al  (1980)  ]  such  transit  time 
effects  are  likely  to  be  important  only  at  very  high  frequencies.  Further, 
the  fairly  complex  spatial  distribution  associated  with  an  actual  dipole 
layer,  as  determined  from  simulations,  suggests  that  a  distinct  transit- 
time  frequency  is  likely  to  be  absent.  Instead  a  spread  in  values  occurs. 

To  summarize;  the  contribution  of  the  drain  noise  component  to  the 
noise  figure  (see  equation  3.3-8): 
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will  show  amplification  within  the  NDM  region  through  ^61  ^  , 
while  the  significant  frequency  dependence  is  through  the  term 
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The  frequency  dependence  is  discussed  at  the  end  of  the  section. 


A  quantitative  evaluation  of  nonlocal  noise  contributions  is  now  con¬ 
sidered  .  For  this  case  intuitive  generalizations  associated  with  the  impedance 
field  method  are  not  as  direct  and  an  alternative  derivation  is  sought.  This 
method  is  first  illustrated  with  an  alternative  calculation  of  the 
nonlocal  contribution  to  6$  following  a  fluctuation.  We  begin  with 
equation  3.4-17  ignoring  the  time  derivative,  a  nonessential  approximation 
that  will  be  relaxed  at  a  later  time  Isee,  e.g..  Richer  (1973)].  Thus 
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where  we  have  again  assumed  Chat  the  dc  mobile  carrier  density  is  uniform. 
Under  open  circuit  conditions  6 id  ■  0  and  a  solution  to  the  resulting 
homogeneous  equation  is  sought.  An  expansion  of  equation  3.4-45  yields 
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Typically  X  is  fo  the  order  of  1Q-^  cm  and  the  contribution  of  the 
second  derivative  term  may  be  safely  Ignored  when  disturbances  are  in 
the  ohmic  region.  Nevertheless,  the  approach  we  are  taking  is  a  new  one 
and  the  discussion  will,  therefore,  concentrate  on  disturbances  within 
this  regime.  The  effects  of  propagation  and  nonlocality  within  the  NDM 
region  will  be  treated  qualitatively,  at  this  time.  We  consider  solutions 
of  the  type  (see  e.g..  Cole  (1968),  section  2.2  on  the  singular  problem). 
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with 
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Rigorously  the  expansion  in  terms  of  X  may  be  justified  only  if  convergence 
is  rapid.  We  are  also  confining  ourselves  to  perturbations  near  one  boundary. 
The  quantities  Sq(x),  S^(x)^,  .  .  .  are  obtained  by  equating  terms  with  equal 
values  of  X,  as  below: 
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Thus 
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and 
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The  problem  is  solved  subject  to  the  condition 
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where  (Xq)  is  presumed  to  be  known.  It  is  then  asserted  that  this 
must  hold  for  any  value  of  X.  Thus 
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Solution  to  equation  3.4-54,  then  yields 


.  0  3.4-i7 


and 


Sc(x) 


A4(Xo) 

ri  -v/*(xo>  -  ♦*] 

ll  -^(X)  -  *g 

3.4-58 


or 


fi40(x)  £  exp  SQ(x) 


3.4-59 


«♦„  <X> 


3.4-60 


41 


Equation  3.4-60  is  the  sane  as  equation  (10)  of  Van  der  Ziel  (1962)  and 
when  used  to  calculate  the  mean  squared  drain  noise,  yields  equation  3.4-32. 
The  first  correction  S^(X)  is  obtained  by  direct  integration  and  is  equal  to 


St(X)  - 


where 


K  ~  f  (+D  "  2  +1  I-*G)  2]  (oex)  -  e(xo))  X 


2L 


3.4-61 


0(X)  0(X  ) 

O 


1  .  0(X)  1  -  0X  ) 


0(X)  -  0(X 


_  tn  f  ._ejX)_.  .  Ll!!V  l) 

,)  l  1  -  «<V  ]) 


0(X)  -  h(X)/H 

(see  equation  3.4-9).  Although  0(X)  >  0(XQ)  for  X  >  XQ,  and  the 
logarithm  term  is  slightly  negative,  (X)  is  negative  and 


X|S1(X)|*  ord  (£) 


3.4-63 


Thus 


and 


d*(X)  *  6^o(X)exp  -  ord 
d*(X)  *  «*0(X)  (l  -  ord 


3.4-64 


3.4-65 


Thus,  for  perturbations  within  the  ohmic  region,  nonlocal  contributions 
to  the  drain  noise  current  are  negligible.  The  situation  when  a  portion  of 
the  FET  is  within  the  NDM  proceeds  in  a  manner  similar  to  that  of  equation 
3.4-50,  except  that  b(X)  is  no  longer  equal  to  unity,  and  a(X)  will  have  a 
form  other  than  that  of  equation  3.4-47,  as  may  be  extracted  from  Bachtold 
(1972).  Further,  X  <  0  and  so  a  disturbance  is  enhanced.  More  generally, 
from  equation  3.4-56  with  given  by  equation  3.4-60  and  S^X)  given  by 
equation  3.4-61,  we  have  the  first  nonlocal  quantitative  correction  to  d^. 

To  calculate  the  short  circuit  gate  noise  it  is  no  longer  assumed 
that  did  “  0.  Instead,  it  is  necessary  to  solve  the  equation 


X  d26*  +  dd±  •  b(X)  +  d#  •  a(X)  -  did  C(X)  3.4-66 
dx2  dx  *L 
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where  a(X)  and  b(X)  are  given  by  equation  3.4-47,  and 

C®  -  (l  -y*<»  -  ♦,  )  3.4-67 

Equation  3.4-66  is  solved  in  a  manner  similar  to  that  of  equation  3.4-45, 
except  that  the  presence  of  the  driving  term  reduces  the  usefulness  of 
the  exponential  expansion.  Instead,  we  let 

64  -  «4C(X)  +  X641(X)  +  X2  642<X)  +  .  .  .  3.4.68 


where 


and 


d«VX)  +  a(X)  6*  (X)  -  did  C(X) 
dx  °  ~ 


dx 


+  a(X)  641(X) 


-d26* 


dx 


3.4-69 


3.4-70 


The  differential  equation  3.4-69  is  the  same  as  that  of  equation  3.4-20. 

If  a  disturbance  in  the  form  of  a  noise  EMF  is  generated  between  X  and  X  +  AX. 

00 

then  a  solution  to  equation  3.4-69  in  the  regions  to  the  left  and  right  of  the 
disturbance,  subject  to  the  short  circuit  conditions 


64(0)  -  64 (L)  -  0 


3.4-71 


yields 


6ijX 


—  “  X  )  (l  -^(X)  -4g) 

Is 

«ld  (X-L)  64o(X)  (l  -  -  4g)) 


0  <  X  <  X  3.4-72 
o 


X  +  AX  <  X  <  L  3.4-73 
o 
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-jiff 


and  within  each  region* 


5*1(X)  “  -e  “  f  a(«)dC  f  dC'e  ^a<C>d4  d  0  < 

J  dc'2 


X  <  X„  3.4-74a 
o 


Z_(x) 

2r  41d 


and 


A  U 

«*1+(X)  “  -e  '  /a(Odt  f  dC’e 


3.4-74b 


fa(t)dr  9 
J  d  d*  (r’) 

'  X  +  AX  <  X <  L  3.4-75a 
•  o 


dC 


,2 


$-r  51d 


3.4-75b 


Across  region  III  the  potential  is  assumed  to  change  by  an  amount  At  (X  ) 
-  o 

Thus 


d$(Xo  +  AX) 

-  d*(Xo)  -  d*o(Xo  +  AX) 

-  dt(X0)  +  X 

fdt1+(X  +  AX)  - 
l  1+  °  3.4-76 

5*1-<V  I 

"  -61d 

1  -yk(X0)  -  H 

•  •  • 

*  81<  \  \ 

3.4-77 

where 

2  »  Z 
A 

ii 

ZI 

3.4-78 

is  independent  of  current.  Thus 

6id  =  -A*(Xo)  | 

i  - 

yMx0)  -  tg  ] 

L  3.4-79 

(i- 

-  *g) 

i 

In  the  X  ■  0  limit,  Vander  Ziel's  (1963)  equation  (11)  is  obtained. 


*The  integration  of  these  equations  is  direct. 
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The  next  step  in  calculating  the  gate  current  noise  involves  comput¬ 
ing  the  fluctuation  in  gate  charge  associated  with  a  fluctuation  in 
potential  A$  at  Xq.  Ve  have  seen  that  such  a  fluctuation  will  produce 
an  alteration  in  the  channel  potential.  If  the  charge  within  the 
depletion  layer  is  changed  by  an  amount 

6Q  -  eNQ  W6h(X)dx  3.4-80 


Then  an  equal  and  opposite  amount  of  charge  will  be  altered  on  the 
gate  contact.  The  net  charge  alteration  is 


3.4-81 


-  -eN  W 
o 


6h(X)dx  -  eNo» 


6h(X)dx  -  eNQW  6h(X)AX 


X  +  AX 
o 


3.4-82 


*  +  AQj  +aqiii  +  aqa 


3.4-83 


liere,  in  Region  I 

<5h  -  Hd» 

2^*  ~  ♦* 

where,  from  equations  3.4-72  and  3.4-74 


3.4-84 


« ♦ 


6ij 


(l  a 


+ 1  Z  (X) 
L  1 


3.4-85 
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For  comparison,  note  equation  3.4-65.  Thus 


AQX  -  -«id 


r  v 

M  e  WH  f 

°  2L  J 


dx*  x* 


y/*  -  *g  (1  V*  "  ♦*> 

+  Vm  j  2i(x)  d*'  | 

°  v/TH7  ) 


3.4-86 


Similarly 


3.4-87 


Thus 


AQ  -  -  6i, 


N  e  WH  J  1  /*  XdX _ 1  f 

-2-j- -  <2  J  0 (X)  (1  -  0<X)>  2  y 


LdX 


8(X)  (1  -  6(X)) 


X  f  Zj.(X)dX  +  *  /• 

2  y  e(x)  2y 


x  +AX 
o 


X  +  AX 
o 


Z_(X)dX  _  1 

-  2 

e(x) 


/ 


X  +AX 

O 


XdX 


0 (X)  (1  -  0(X) 


-e  N  WH  A*  AX 
o  _  o 

— T~  e(xo) 


3.4-88 
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where  8(Xq)  Is  given  by  equation  3.4-62.  Integration  of  the  first  two 
terms  Is  direct.  Using  equation  3.4-19  along  with 


3/ 

X  -  -  ^  |«(X)  -  |  j  (♦«)  -  4g)  2  -  <-9g>  "I  3.4-89 


3/,) 


we  obtain 


AO  »  -6i .  N  eLUH  j®  "  0(Xo*  +  X  Z„(X_)( 

d  °  |“5T"  “  T  ) 


3.4-90 


where 

Z(V  '  E 

f Zx(X)dX 

+  /  ZTT  (X)dXl 

J  1  e(x) 

Xo+& 

3.4-91 

fo  9<» 

> 

and 

L 

l 

L 

ye 

/ 

5-i  f 

e(x)dx 

3.4-92 

o 


The  first  part  of  equation  3.4-90  we  recognize  as  the  same  as  equation  (16) 
of  Van  der  Ziel  (1963). 

We  are  now  in  a  position  to  compute  the  mean  square  noise  currents. 
Recognizing  that  A$(Xq)  is  caused  by  thermal  noise,  its  mean  squared 
value  is: 


d*(XQ) 


4kT  Af  AX 


dl 


dtp 


d(X0)  L 


3.4-93 


or  in  normalized  units,  as 


A^(X(J)2-  4kT  Af  •  AX 
2. 


G  ♦„  Ld  Ao/ d$ 
on  A 

Hence,  for  the  drain  current  noise  (from  equation  3.4-79): 


3.4-94 


fiij 


4kTAf • AX 


>VH2dldo*°,/d* 


.  /  1  -  /♦«„>  -  ♦«  \ 
\  1  +  XZA(1  -^*(Xo)  -  *g)  / 


and 


<K>  ‘  / 


6id  dx 


Noting  that 


dMxJ  / - 

-Sr*-  -  1  -  -  ♦ 


g 


equation  3.4-95  integrates  to 


3.4-95 


3.4-96 


3.4-97 


<«l^>  *  «IM_  |«W-i  /  (l  ZA(X>  to}  3.4-98 


oTH 

Or  in  ordinary  units 


3 

4kIG°  /  f1  -V*<v  -  ♦«)  VV*] 


3.4-99 


whereP($,$)  is  given  by  equation  3.4-32.  We  note  the  similarity  of 
D  G 

equation  3.4-99  to  the  more  intuitive  argument  of  equation  3.4-41. 
For  the  mean  squared  fluctuating  gate  charge 


AQ  -  (NQeLWH)‘ 


e  -  e(xo)  +  xz(xo) 


6ij 


3.4-100 


48 


-  <NoeLWH)  I  S  -  9<Xo)  +  XZ(XD 


do 


1  -^(xo)  -  *8 


'1  - 
1  + 


AZa(1  -yj(Xo) 


) )  •  4kTAf  AX 

-f  ir  *  2  dl^<x 

)  lg0*h  -Sr 

of 

2 

-  $g| 


The  mean  squared  charge  fluctuation  along  the  entire  gate  is 

L 


<*»  ■  / 


77T2  d* 

AQ 


Defining 


Csg  - 


N  eLWH 
o 


WV 


r  ^1  -y»(X)  -  *8  j  d* 

A 


*d  _  /  _ \2 

h2(«d.«g)  -  y*  ^(x)  -  *g  y*(X)  -  *g  j 


o 

♦a 


h3  <4d,4  > 


/ 


(♦(X)  -  *g)  -y*(X)  -  *g~J  d4 


[see  Vander  Ziel  (1963)]  with  (see  equation  3.4-32) 


fl“  Ldo 

the  uncorrected  portion  of  3.4-102  integrates  to 

<*v  -  •  c. J  •  <jvy 

wy3 

where 

g3(4d*V  “  ®2  W*g}  "  2®h2(*d’V  +  h3(*d’V 


3.4-101 


3.4-102 


3.4-103 


3.4-104 


3.4-105 


3.4-106 


2. 4-107 


3.4-108a 


3.4-108b 


The  total,  including  corrected  contributions,  is 
2 


_  4kT  Csg  •  )  £. 


Afl/2»  G 


-  L 

/* d*(i-e)2(e 

l<*d»0  L 

8  -  L  o  i 


-e)z(xo) 

2 

do 


-  2Za(1-8)3  (6-0) 2 


do 


The  mean  square  gate  current  is 


g' 

For  the  correlation  coefficient 


<«i  >2  -  n»2 


or 


and 


AQ«id  -  NQeLWH 


e-e(x  > 


*do  2L 


'SL  +  _L  20O  > 

2L  °  J 


AQ6Id  -  AQ6id  GJh 


AQ6I , 


VLWHGo»H  ’4KJAf  AX 

Go*h2L  ido/ 


(*-e(V  +  »*<V)  1 1  -  /^xo)  "  * 

— n- 


1  +  -  ViTxIJ  -  4>g 


Thus, 


la 

<fo6Id>  /" 

dn/2n  "  y 


3.4-109 


3.4-110 


3.4-111 


3.4-112 


3.4-113 


o 


AQ«Id  dx 


3.4-114 


*  (i  /  ° 


where 


3.4-115 


-l  fzt  l  -  Vmo) $g)3  (e-e(x»dJ 


g2(4d.O  "  ®hi  **d»*g^  ~  h2  ^d**g* 


3.4-116 


Mote  the  first  part  of  equation  3.4-115  Is  found  In  Vander  Ziel 
(1969) .  If 


61  -  -JO  AQ 

8 


3.4-117 


<«.  «d>  •  -jn  <AQ6Id> 


3.4.5  Conclusions 

Equations  3.4-99,  3.4-109  and  3.4-115  represent  the  first 
significant  analytical  generalization  of  FET  noise  theory  to  account 
for  nonlocal  effects  of  space  charge  within  the  conduction  channel  of 
the  FET.  The  application  was  only  to  ohmic  conduction,  but  an  apparent 
generalization  to  transport  within  the  NDM  region  is  clear.  The  impor¬ 
tance  of  this  calculation  lies  in  the  fact  that  the  current  sources 
depend  significantly  on  the  space  charge  distribution  within  the  FET, 
and  that  disturbances  within  the  NDM  regime  will  be  amplified. 

The  effect  of  coupling  this  to  the  external  circuit  is  unclear. 
Consider  for  example  equation  3.4-44  which  represents  the  effect  of  drain 
circuit  and  noise  source  on  the  noise  figure  of  the  FET.  At  low  values 
of  drain  bias,  figure  3-9  indicates  that  at  all  but  dc  conditions,  the 
imaginary  part  of  Y21  is  the  dominant  contribution.  Similarly,  for 
the  imaginary  part  of  Y^.  But,  in  this  case  there  is  no  enhanced 


noise  due  to  gain  from  within  the  NDM  region.  Nevertheless,  under  dc 
conditions,  the  situation  appears  as 


3.4-118 


Equation  3.4-118  prevails  at  high  bias  levels  when  domains  are  present 

and  introduces  two  competing  contributions.  On  the  one  hand,  increased 

2 

noise  is  expected  from  <61^  >,  but  this  is  partially  compensated  by  an 

increase  in  the  transconductance.  Thus,  the  net  effect  may  be  an  apparent 
decrease  in  noise  due  to  the  drain  current  contribution. 

The  frequency  dependence  at  bias  levels  sufficiently  high  to  generate 
domain  is  more  direct.  Here,  Y^  is  dominated  by  capacitive  contributions 
while  Y2^  is  still  primarily  real.  Thus,  the  drain  component  of  noise  is 
approximately 

<«£>  1  »  ■*  T.2  Cn  3.4-119 

gm 


and  shows  an  increase  with  increasing  frequency,  as  discussed  experimentally 

[Liechti  (1976)  ].  As  far  as  the  frequency  dependence  of  the  gate  noise 

current,  there  is  no  contribution  from  the  circuit  as  long  as  Ys  is  real. 

2 

All  frequency  dependence  is  contained  in  <  61^  >.  With  regard  to  the 
correlation  coefficient,  its  frequency  dependence  is  on  the  square  root 
of  the  drain  component. 
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Figure  3-1.  Two  dimensional  device  circuit  configuration  for  the  small  signal  calcula¬ 
tions.  The  semi-conductor  studied  was  gallium  arsenide  with  a  nominal  doping  Nq  «•' 
lO^/CM?.  The  low  doping  was  chosen  to  reduce  computational  time.  Relevant  dimensions 
are  1^  ■  1.2pm,  LgQ  *  10pm  and  H  -  1.95'jm. 


Figure  3-2.  Steady  state  drain  current  versus  drain  potential  for  two  values  of 
gate  bias.  Here  Iq  ■  NQeV  A  and  \l>0  -  Fpl,  where  V_  is  the  peak  gallium  arsenide 
carrier  velocity,  and  F_  is  the  threshold  field  for  negative  differential  mobility. 
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CHANNEL  LENGTH 


3-4.  Carrier  density  showing  a  broad,  almost  symmetric,  carrier  density.  Por  this  calculation 
.6^0,  and  “  0.44iq.  The  current  density  for  this  case  is  In»0.0. 


Figure  3-5.  Representation  of  the  small  signal  dlstrubance.  A  square  wave  pulse 
Is  applied  to  the  drain  contact,  altering  the  drain  contact  boundary  condition: 

4*d  “xf’D  +  fi'f'D-  There  Is  a  subsequent  change  in  the  drain  and  gate  currents, 
dip  and  dig. 


- 

Bras 

'  - '» 

.  ? 

Figure  3-6  -  General  form  of  the  FET  Intrinsic  circuit  model  with  reference  to 
equations  3.2-11  and  3.2-12;  *  ?22  +  *12*  *11  ■  *n  +  *12* 

*21  "  -*12  *«d  *ml  -  Y21  “  *12.  See  also  Cobbold  (1970),  figure  5.3. 
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Figure  3-7  -  First  order  equivalent  circuit  for  the  small 
signal  FET  calculations. 


Figure  3-8  -  Schematic  representation  of  the  gate-drain  admittance  in  the  presence 
of  a  high  field  domain.  Gg«j  represents  a  small  signal  conductance 
which  exhibits  a  frequency  dependent  region  of  negative  differential 
conductivity  when  a  high  field  domain  is  present. 
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PARAMETERS 


3.9b  -  Frequency  dependent  admittance  Yu  and  Y21  obtained  by  Fourier 
analysis  of  gate  and  drain  currents,  when  the  gate  is  subjected 

to  a  snnarp  wavp  ntilap. 


ELECTRON  DENSITY 


state  electron  density  for  this  bias  level 
discussion  of  projection  characteristics. 


PARAMETERS 


FREQUENCY 


FREQUENCY  FREQUENCY 


PARAMETERS 


FREQUENCY 


PARAMETERS 


Figure  3-13a,  b  and  c  -  Aa  in  figure  3-9,  but  with  Fjjo  -  0.1Fo  and  TGo  ■  -0.4T, 


PARAMETERS  Fig.  3-140 


FREQUENCY  FREQUENCY 


Fig.  3-l4b. 


CAPACITANCE 


(^00  "  *oo)/*0 


Figure  3-15  -  Gate  to  drain  capacitance  versus  (Tug  -  Vqo^^o'  Circular  and 
triangular  results  are  for  a  fixed  gate  bias.  Square  and 
inverted  triangular  data  are  for  a  fixed  drain  bias.  For  this 
calculation  Cq  *  3.25  x  10^  cHW/L. 
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Figure  3-18  -  Capacitance  (see  equation  3.2-41)  versus  (^dq  -  '•'go^^o- 
For  this  calculation  Cq  =  eHW/L.  Note  »  Cgs  +  Cgd- 


■1  •  }  «  *  •  It 


v^/v  — — 

Figure  3-19  -  Experimental  measurements  of  the  gate-drain  capacitance 
From  Enjslenann,  et  al  (1977). 


Figure  3-20  -  Re  Yj^  (in  multiples  of  GQ)  versus  (Yjjq  -  ’•'go^^o*  There  is 
saturation  at  high  bias  levels.  The  slight  decrease  in 
transconductance  at  high  bias  levels  for  YGq  ■  -0.lYo  may  be 
a  numerical  artifact.  It  is  being  studied. 


Figure  3-21  -  Experimental  measurement  of  transconductance  versus  drain  bias 
From  Englemann,  et  al  (1977). 


CUTOFF  FREQUENCY  (GHZ) 


Figure  3-23  -  Experimental  measurements  of  the  cutoff  frequency  versus 

From  Englemann,  et  al  (1977). 
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Figure  3-24  -  Circuit  representation  of  an  FET  with  noise. 


Figure  3-26  -  Circuit  representation  of  a  noisy  FET  as  a  noiseless  network 
with  a  voltage  and  current  source.  [See  Talpey  (1959)]. 


Figure  3-27  -  As  in  figure  3-26,  but  with  signal  noise  source.  [From  Talpey  (1959)]. 
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Figure  3-28  -  Structure  of  FET  and  definitions  of  variables  used  in  gradual 
channel  approximation  steady. 
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4.  Thermal  Variations  in  Semiconductor  Devices 


4.1  Introduction 

When  a  material  such  as  gallium  arsenide  or  indium  phosphide  is  operating 
as  a  Gunn  diode,  where  the  operation  is  dependent  upon  transit  of  a  high 
field  domain,  the  resulting  frequency  is  dependent  in  a  significant  way  on  the 
saturated  drift  velocity  of  the  semiconductor.  The  saturated  drift  velocity  of 
the  carriers  is  significantly  dependent  upon  the  high  field  electron-phonon 
scattering  rates  and  in  a  marginal  way  upon  scattering  by  ionized  Impurities. 

At  high  fields  the  saturated  drift  velocity  decreases  with  increasing  temperature, 
as  seen  in  figure  4-1.  In  addition  to  changes  in  the  high  field  temperature, 
dependent  properties  of  the  velocity-field  curve  there  is  also  a  decrease  in  the 
low  field  mobility  with  increasing  temperature.  Thus,  it  may  be  expected  that 
increases  in  the  operating  temperature  will  degrade  device  performance. 

If  we  move  away  from  the  Gunn  device  and  toward  a  compound  semiconductor 
IMP ATT  it  becomes  necessary  to  include  hole  as  well  as  electron  transport. 
Furthermore,  temperature  variations  in  a  direction  normal  to  the  principle 
direction  of  current  flow  are  known  to  cause  measurable  changes  in  the 
performance  of  IMPATTS  [Haitz  (1968)].  Thus,  two-dimensional  simulations  are 
necessary.  The  two-dimensional  nature  of  the  problem  becomes  even  more 
dramatic  when  it  is  realized  that  gallium  arsenide  FETs  and  three- terminal 
transferred  electron  logic  devices  are  extremely  sensitive  to  temperature  varia¬ 
tions  [Grubin,  et  al  (1980,  1982)].  Under  the  present  ARO  contract,  two-dimensional 
transient  variations  for  both  hole  and  electron  transport  were  formulated.  However, 
numerical  calculations  we  performed  for  only  one  space  dimension,  and  it  is 
within  this  context  that  our  results  are  presented. 

The  temperature  dependent  results  discussed  below  must  be  regarded  as 
initial  effects  on  this  problem.  Additional  calculations  are  still  necessary 
before  preparation  for  publication. 


4.2  Formulation  and  Analytical  Results 


Thermal  effects  are  examined  by  generalizing  equation  2.23  to  include 
temperature  effects: 

J(T)  -  NeV(E,^  -  e  ^(E,^)  3N^*-T)-  +  D(E,£*)N(X,T)  fjH  4.2-1 


+  c  ss. 

e  3T 


C)  -  e.  D 
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The  third  tern  on  the  right  side  of  Eq.  4.2-1  represents  temperature  gradient 
contributions  to  the  current  density.  It  is  included  for  formal  reasons  only, 
and  by  analogy  to  the  situation  in  which  classical  statistics  apply  and  the 
scattering  parameter  is  Independent  of  position.  Its  contribution  will, 
however,  be  ignored.  The  quantities  considered  below  are  the  field-dependent 
velocity  and  diffusion  coefficients  which  now  Include  a  temperature  dependence. 
The  temperature  dependence  of  the  velocity-electric  field  curve  has  been 
calculated  by  Ruch,  et  al  (1970)  and  shown  to  be  in  good  agreement  with  the 
experimental  measurements  of  Ruch,  et  al  (1968).  A  fit  to  the  velocity-field 


curve  of  Ruch,  et  al  (1970)  was  given  by  Freeman,  et  al  (1972) 

@)‘]l 


4.2-2 


and  is  said  to  be  valid  for  any  temperature  between  300  and  600°K.  The  fit  for 
300  and  500*K  is  shown  in  Fig.  4-1.  In  Eq.  4.2-2,  Eq  is  4  k  V/cm,  yielding 
a  threshold  electric  field  of  approximately  3.5  k  V/cm.  Freeman,  et  al  (1972) 
also  argued  that  the  temperature  dependence  of  the  diffusion  curve  between  300 
and  600°K  could  be  represented  by  multiplying  the  room  temperature  values  by  the 
factor  300 °K/  Analytical  representations  of  the  velocity  field  curve,  as 

given  by  equation  4.2-2  are  the  type  used  in  the  simulation. 

To  simulate  the  temperature  dependence  of  the  device,  equation  4.2-1  was 
solved  simultaneously  with  the  equation  for  heat  conduction 

CPP  Uf)  ffU  J(T)E(X,T),  4.2-3 


where  p  is  the  mass  density  of  the  material  and  C  its  specific  heat  [see  Carslaw 
and  Jaeger  (1959)].  In  studying  this  problem,  two  time  constants  emerge.  The 
thermal  time  constant  [see  e.g.,  Shaw,  et  al  [1979]  equation  7-30) 


t.  -  4Cp*  2/r2<  «  (1.5  x  *  2/cm2)s,  4.2-4 

tn  n  n 

and  a  characteristic  electric  time  constant.  For  1  <*  lOym  r  ,  -  1.6  x  10  ^  sec. 

n  tn 

In  the  one-dimensional  gallium  arsenide  calculations,  the  doping  profile  were 
chosen  to  yield  a  localized  high  field  domain  which  was  potentially  capable  of 
sustaining  a  time-dependent  oscillation.  At  sufficiently  high  bias  an  oscillation 
did  result.  The  time  constants  associated  with  this  oscillation  were  of  the 


85 


order  of  £  /V,  and  for  the  device  studied  £  -  Sum  leading  to 
n  n 

rel  *  £q/V  *  5  x  ltf11  sec  *  3  x  1<T5  t 


4.2-5 


The  results  are  expected  to  lead  to  a  slow  Increase  in  temperature  from  the 
source  to  the  drain.  However,  in  view  of  the  disparate  values  of  the  thermal 
and  electrical  time  constants,  we  chose  not  to  resolve  the  long  thermal  time 
to  equilibrium.  Instead,  we  concentrated  on  the  steady  state  time-independent 
solutions.  Before  discussing  these  solutions,  it  is  useful  to  make  reference 
to  an  earlier  analytical  calculation,  first  discussed  by  Knight  (1967)  and 
summarized  by  Shaw,  et  al  (1979). 

We  consider  the  device  structure  shown  in  figure  4-2  and  seek  time- 
independent  solutions  to  equation  4.2-3  for  the  situation  when  the  quantity 

4>  i  I 

of  heat  generated  at  the  n  n  Interface  and  within  the  n  region  is  negligible, 
are  [Shaw,  et  al  (1979)] 


-  rn/n-H-  exp  ^|an2-x2)},  0  <  X  <  ln. 


4.2-5 


where  ^n/n^  is  the  temperature  at  the  n/n  interface  and  £q  is  the  thickness 
of  the  n  region. 


^•(X)  -  ^  exp  |uEtn/120) (£n  +  £n++  -  X)]  , 


4.2-6 


£  <  X  <  £  ++  +  £ 

n  n  n 


where  £  ++  is  the  thickness  of  the  n  layer  and  v  is  the  temperature  at 
n  ms 

the  metal/semlconductor  interface.  From  equation  4.2-5  and  4.2-6 


4.2-7 


and  the  temperature  distribution  within  the  n  layer  is 

T(X)  -  r,  expJE£n  (^ko-  +  3^)]  [***»  fw)]  ’ 


0  <  X  <  £ 


We  are  interested  in  the  temperature  difference  ^(0)  -  -  if 


4.2-8 
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Thus 


tf  -  f 

w  v  ms 

% 

If  1  +f  <  <  i 
n  n 

“  Cis 

For  the  parameters  of  Table  2,  and  a  power  density  associated  with  the  maximum 

field  and  current  prior  to  an  Instability,  JE4q  ”  5.75  x  10^  W/cm^, 

and 

5  ^  101]  4 . 2-11 

For  m  300®,  A^** “  3“K  for  a  field  of  90KV/cm 

IDS 

lC  =  C  [0.11]  4.2-12 

m 

and  again  for  *  300°,  A  =  34®K 

Thus,  for  very  short  devices  large  temperature  gradients  are  not  expected  until 
high  fields  are  reached.  These  same  general  conclusions  emerge  from  the 
numerical  calculation.  The  results,  including  the  doping  profile  are  shown  in 
figures  4-3  through  4-6.  With  the  exception  of  one  calculation,  all  were 
performed  for  a  thermal  conductivity  of  X  ■  0.5  W/cm°K.  The  temperature 
dependent  boundary  conditions  were:  @  X  ■  0,  2^  »  300  “K,  @  X  «  L,  3  /3X  *  0. 

4.3  Numerical  Calculations  and  Conclusions 

Figure  4-3  shows  the  carrier  density,  electric  field  profile,  and  temperature 
variation  for  a  5pm  long  gallium  arsenide  element  with  a  centrally  placed  notch 
of  magnitude  An  •  .6Nq.  The  potential  drop  across  this  element  is  1  volt. 

Note  the  other  negligible  temperature  variation  across  the  device.  The  current 

level  for  these  calculations  yielded  an  average  velocity  <  v  >  »  J/N  e  * 

7  0 

1.46  x  10  cm/sec.  Further  increases  in  voltage  generated  transit  time 

oscillations. 

Figure  4-4  displays  a  similar  calculation.  This  time,  however,  with  a 
doping  notch  of  An  ■  0.1Nq.  For  this  case,  the  field  variation  is  more 


extreme  than  in  figure  4-3,  but  the  current  level  is  below  that  of  the  former; 
a  result  consistent  with  the  cathode  boundary  field  model  [See,  e.g.,  Shaw,  et  al 
(1979)].  Note  that  for  both  figure  4-3  and  4-4  the  temperature  variation  is 
negligible,  as  expected  from  equation  4.2-11.  With  regard  to  figure  4-4,  we 
point  out  that  the  structure  in  N(X) ,  downstream  from  the  notch  is  a  consequence 
of  the  region  of  negative  differential  mobility. 

Large  temperature  variations,  for  a  sample  of  fixed  length,  are  a 
consequence  of  either  a  low  thermal  conductivity  and/or  high  fields.  We 
illustrate  both  cases.  To  illustrate  the  low  thermal  conductivity  case,  we  use 
gniHiim  arsenide  parameters  with  a  thermal  conductivity  equal  to  0.05  W/em  *K. 

The  result  is  shown  in  figure  4-5.  While  the  current  density  is  virtually 
unchanged,  the  temperature  variation  has  increased  by  an  order  of  magnitude.  . 

The  situation  in  which  an  average  field  is  increased  to  90  kv/cm,  and  the 
thermal  conductivity  taken  as  0.5  W/cm  °K,  is  displayed  in  figure  4-6.  While 
structure  in  carrier  density  and  field  is  seen  everywhere  throughout  the  device, 
the  field  profile  is  relatively  uniform.  For  this  calculation,  the  origin  of 
all  field  nonuniformities  is  the  notch,  and  we  see  that  the  temperature  dif¬ 
ference  between  cathode  and  anode  is  approximately  30 “K.  This  was  predicted 
analytically  in  equation  4.2-12. 

While  additional  calculations  are  necessary,  particularly  with  regard  to 
ellucldatlng  the  role  of  field  nonuniformities  on  the  temperature  profile,  the 
results  of  both  the  analytical  and  numerical  study  Indicate  that  for  short 
micron  scale  devices  large  thermal  gradients  are  not  likely  to  be  prominent 
until  high  fields  are  present  across  the  device  structure. 


TABLE  2 


V  -  2.25  x  107  cm/sec 
P 

E  -  3.2  kV/cm 
P 

-4 

t  ■  J  x  10  cm 

N  ■  1013/cm3 
o 

-4  2 

area  *  1.5  x  10  cm 


GaAs  Parameters 


R  -  3.030 
o 

C  «  2.94  x  10"13f 
o 


J  -  N  eV 
p  op 


-  3.6  x  103  A/cm2 


El*  1.6  volts 
P 


Figure  4- 


Figure  4-: 


*  . 


ELECTRIC  FIELD  (W/CM-1l 


1  -  Analytical  fit  ( - )  of  the  drift  velocity  electric  field 

relation  for  electrons  in  GaAs.  (From  Freeman  et  al  (1972). 
For  the  data  ( — ) ,  see  Ruch  et  al  1968 . 
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-  Device  structure  for  doing  temperature-dependent  calculations. 
(From  Knight  (1967). 
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CARRIER  DENSITY  (N/(^) 


Figure  4-3(a)  -  Normalized  carrier  density,  N/lO^Van3.  and 
normalized  background  doping  No(x)/10^^/cm3 
for  a  five  micron  long  gallium  arsenide  element 
subjected  to  an  average  field  of  2kv/cm.  The 
thermal  conductivity  is  ic  ■  0.5W/cra*K. 
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figure  4-3 (c)  -  Fractional  temperature  variation  across  device. 
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□STANCE  (microns) 


Figure  4-4  a,b  and  c  -  As  in  Figure  4-3,  but  for  a  different 

notch  depth. 


HELD  (HV/  cm) 


Figure  4-5  a,b  and  c  -  Aa  in  Figure  4-4,  but  for  a  reduced 

thermal  conductivity,  k  ■  0.05W/cm*K 


DISTANCE  (microns) 


Figure  4-6 


a,b  and  c  -  As  in  Figure  4-4,  but  for  an  average 
field  of  90kv/cm. 
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5.  Drift  and  Diffusion  Description  of  Injection  Field 
Effect  Transistors 


5.1  Introduction 

A  recent  study  by  Fauquemberque,  et  al  (1983)  on  a  submicron  gate 
injection-type  field  effect  transistor  suggests  some  rather  interesting 
electrical  characteristics.  In  particular,  for  a  gate  length  of  O.lSvm,  a  gain- 
bandwidth  product  of  400  GHZ  was  predicted.  The  predicted  behavior  of  the  device 
was  based  on  a  Monte  Carlo  solution,  and  short  channel  velocity  overshoot 
effects  were  included.  The  inmediate  question  is:  To  what  extent  are  submicron 
dimensions  important  for  this  pehnomena.  A  study  to  answer  this  question 
was  recently  initiated.  The  initial  study  was  for  a  micron  scale  device. 

The  calculations  were  performed  for  a  device  whose  basic  structure  is 

shown  in  figure  3-1.  However,  the  semiconductor  material  used  for  this 

16  3 

calculation  was  InP  with  a  nominal  doping  level  of  10  /cm  ,  H  «  0.64  micron 
and  L  ■  3.0  microns. 

5.2  numerical  Results 

Two  types  of  calculations  were  performed.  In  one  calculation  was  uniform 
and  at  sufficiently  high  bias  the  device  sustained  current  oscillations  in  the 
form  of  propagating  high  field  domains .  In  the  second  calculation,  the  back- 
ground  was  reduced  to  10  /cm  over  a  distance  of  1.5vm.  Ho  instabilities  were 
observed  in  the  latter  case  and,  of  course,  the  current  levels  were  significantly 
reduced.  The  space  charge  distribution  for  the  uniformly  doped  material  was 
qualitatively  similar  to  that  obtained  for  the  gallium  arsenide  structure 
discussed  in  section  3.  The  distribution  for  the  injection  FET  is  qualitatively 
different,  and  is  displayed  in  figure  5-1,  where  we  see  strong  injection. 

The  first  collection  of  calculations  for  the  injection  FET  provide  a 
series  of  current  voltage  characteristics,  two  branches  of  which  are  shown  in 
figure  5-2.  The  initial  results  find  a  dramatic  Improvement  in  the  trans¬ 
conductance  for  the  injection  FET.  For  example;  For  the  uniformly  doped  FET 
and  a  drain  potential  of  approximately  0.5  volts  the  transconductance  near 
▼G  -  -0.1  was 

5.2-1 


gm  “  .5 

o 
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where 


N  epHW 

Go  -  ■  •  L-  -  5.2-2 

For  the  Injection  FET  Gq  Is  reduced  by  almost  two  orders  of  magnitude: 

Golnj  «  G  /50  5.2-3 

o 

while  the  transconductance  is  reduced  by  significantly  less 

gn  *  0.14  Gq  5.2-4 

The  results  of  these  studies  are  currently  being  evaluated,  but  several  points 
are  noted: 

For  a  routinely  configured  FET,  the  change  in  drain  current  due  to  changes 
in  gate  bias  is  largely  a  consequence  of  a  modulation  in  the  cross-sectional 
area  of  the  conducting  channel.  When  domains  form  as  discussed  in  section  3, 
modulation  is  accompanied  by  space  charge  injection  into  the  depletion  region. 
For  the  injection  FET  the  space  charge  density  at  moderate  value  of  gate  bias 
levels  exceed*  the  background.  Modulation  is  then  accompanied  by  both  a 
variation  in  channel  height  as  well  as  an  alternation  in  the  net  mobile  carrier 
density.  The  behavior  of  the  injection  FET  appears  operationally  different  from 
the  classical  FET,  and  its  device  potential  should  be  explored. 


Figure  5.1a  -  Spatial  distribution  of  carrier  density  N(X)  at  the 
bottom  of  the  channel  of  the  injection  FET.  Also  shown  is  the  back¬ 
ground  density  Nq(X)  for  the  device. 


'K'f'o 


Figure  5.1b  -  Potential  distribution  along  the  bottom  of  the  channel. 

Most  of  the  potential  drop  is  across  the  Injection  region,  as  shown  by 

the  solid  line  whose  axis  is  toward  the  right  hand  side  of  the  figure. 

There  is  structure  in  the  potential  near  the  source.  This  is  shown  by 

the  dashed  curve.  For  this  calculation  T  ■  3  volts. 

o 
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